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Abstract 

In this paper we consider random block matrices, which generalize the general beta en- 
sembles, which were recently investigated by Dumitriu and Edelmann (2002, 2005). We 
demonstrate that the eigenvalues of these random matrices can be uniformly approximated 
by roots of matrix orthogonal polynomials which were investigated independently from the 
random matrix literature. As a consequence we derive the asymptotic spectral distribution 
of these matrices. The limit distribution has a density, which can be represented as the 
trace of an integral of densities of matrix measures corresponding to the Chebyshev matrix 
polynomials of the first kind. Our results establish a new relation between the theory of 
random block matrices and the field of matrix orthogonal polynomials, which have not been 
explored so far in the literature. 
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1 Introduction 



The classical orthogonal, unitary and symplectic ensembles of random matrices have been studied 
extensively in the literature, and there are numerous directions to generalize these investigations 
to the case of random block matrices. Several authors have studied the limiting spectra of such 
matrices under different conditions on the blocks [see e.g. Girko (2000), Oraby (2006 a,b) among 
others] . 

In the present paper we study a class of random block matrices which generalize the tridiagonal 
matrices corresponding to the classical /3-Hermite ensemble defined by its density 



l<i<j<n 



En 



(1.1) 



where /5 > is a given parameter and C/? a normalizing constant. It is known for a long time that 
for /3 = 1, 2,4 this density corresponds to the Gaussian ensembles [see Dyson (1962)], which have 
been studied extensively in the literature on random matrices [see Mehta (1967)]. It was recently 
shown by Dumitriu and Edelmann (2002) that for any f3 > the joint density of the eigenvalues 
Ai < ■ ■ ■ < A„ of the n X n symmetric matrix 

N2 



75^{n-l)/3 



n-l 



V 



Nn J 



;i.2) 



is given by fll.ip where X^, . . . , X(^n-i)f3, Ni, . . . ,Nn are independent random variables with Xj^ ~ 
X^f^ and Nj ~ A/'(0, 1). Here A'/^ and A/'(0, 1) denote a chi -square (with jjS degrees of freedom) 
and a standard normal distributed random variable, respectively. In the present paper we consider 
random block matrices of the form 



^ip) q{p) ^(p) 
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:i.3) 



^^-1 / 



where n = mp G N, m G N, p G N and the p x p blocks Af'^ and B^^^ are defined by 
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ip+l 
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7i(n— jp— 1) ^72(n-— ip— 2) 
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71 (n—ip—2) 



X. 
X. 



'yp-i{n-{i+l)p+l) 
7p-2('T.-{«+l)p+l) 



-^7p_2(n-(i+l)p+2) 
y -^7p_i(n-(i+l)p+l) 



X. 



7l(n-(i+l)p+l) 



V2N( 



(i+i)p-i 



X. 



7i(n-(i+l)p+l) 



-^72(n-(j+l)p+l) -^7i(n-(i+l)p+l) 



(1.4) 



(i+i)p / 



and 
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(1.5) 



respectively, 71, . . . ,7p > are given constants and all random variables are independent. Note 
that in the case p = 1 we have Gn^ = Gn^ and that for p > 1 the elements g^^j^ of the matrix Gn^ 
are standard normally distributed if i = j and Af^-distributed else, where the degrees of freedom 
depend on the position of the element in the matrix Gn^ . In the present paper we investigate the 
limiting spectral behaviour of matrices of the form ( 11. 3p . In particular it is demonstrated that the 
eigenvalues of the matrix Gn^ are closely related to roots of matrix orthogonal polynomials, which 
have been studied independently from the theory of random matrices [see e.g. Duran (1996, 1999), 
Duran and Daneri-Vias (2001), Duran and Lopez- Rodriguez (1996), Duran, Lopez- Rodriguez and 
Saff (1999), Sinap and van Assche (1996) or Zygmunt (2002) among others]. To our knowledge 
this is the first paper, which relates the field of matrix orthogonal polynomials to the theory of 
random matrices. 

In Section 2 we review some basic facts on matrix orthogonal polynomials. We prove a new result 
on the asymptotic behaviour of the roots of such polynomials which is of own interest and the 
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basis for the investigation of the hmiting spectral properties of the matrix Gn ■ In Section 3 
we provide a strong uniform approximation of the random eigenvalues of the matrix Gn^ by the 
deterministic roots of matrix valued orthogonal polynomials, and these results are applied to 
study the asymptotic behaviour of the spectrum of the matrix G^^. In particular, we derive a 
new class of limiting spectral distributions which generalize the classical Wigner semi circle law 
in the one-dimensional case. Roughly speaking, the limit distribution has a density, which can 
be represented as the trace of an integral of densities of matrix measures corresponding to the 
Chebyshev matrix polynomials of the first kind. Finally some examples are presented in Section 
4 which illustrate the theoretical results. 

2 Matrix orthogonal polynomials 

In the following discussion we consider for each A; e N a sequence of p x p matrix polynomials 
{Rn,kix))n>o, which are defined recursively by 



with initial conditions R^i^k{x) — 0, Ro^ki^) — Ip, where B^k £ denote symmetric and 

Ai,k e K^^^ denote non-singular matrices. Here and throughout this paper Ip denotes the p x p 
identity matrix and a p x p matrix with all entries equal to 0. A matrix measure E is a p x p 
matrix of (real) Borel measures such that for each Borel set A C M the matrix G M^'^^ is 

nonnegative definite. It follows from Sinap and van Assche (1996) that there exists a positive 
definite matrix measure e W^^^ such that the polynomials (-Rn,fe(2;))„>o orthonormal with 
respect to the inner product induced by the measure dEfe(x), that is 



The roots of the matrix polynomial Rn,k{x) are defined by the roots of the polynomial (of degree 
np) 



and it can be shown that the orthonormal matrix polynomial Rn,k{x) has precisely np real roots, 
where each root has at most multiplicity p. Throughout this paper let Xn,k,j {j — ^, ■ ■ ■ i'itt) denote 



xRn,k{x) = A„+i^kRn+l,k{^) + Bn,kRn,k{x) + A^,^Rn-l,k{x) , 



(2.1) 




(2.2) 



deti?„,fc(x). 
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the different roots of the matrix orthogonal polynomial Rn,k{x) with multiplicities ij, and consider 
the empirical distribution of the roots defined by 

^ m 

^"'^ •= ~ Zl^J'^w.^ ' ^' ^ ^ 1' (2-3) 

where denotes the Dirac measure at point z G M. In the following we are interested in the 
asymptotic properties of this measure if n, /c ^ oo. For these investigations we consider sequences 
(^j)jeN) (^j)jeN of positive integers such that ^ ^ m for some u G [0, oo), as j ^ oo and denote 
the corresponding limit as lim„/fc^„ (if it exists). Our main result in this section establishes the 
weak convergence of the sequence of measures 6n,k under certain conditions on the matrices An,k 
and Bn,k if n/k —>■ u (as n oo, k —>■ oo) in the above sense. In the following sections we will 
approximate the eigenvalues of the random block matrix Gn^ by the roots of a specific sequence of 
matrix orthogonal polynomials and use this result to derive the asymptotic spectral distribution 
of the random block matrix Gn\ 

The theory of matrix orthogonal polynomials is substantially richer than the corresponding theory 
for the one-dimensional case. Even the case, where all coefficients in the recurrence relation (12. ip 
are constant, has not been studied in full detail. We refer the interested reader Aptekarev and 
Nikishin (1983), Geronimo (1982), Dette and Studden (2002) and the references cited in the 
introduction among many others. Before we state our main result regarding the asymptotic zero 
distribution of the measure 6n,k we mention some facts which are required for the formulation of 
the following theorem. The matrix Chebyshev polynomials of the first kind are defined recursively 

by 

tT^'^it) = AT^'^{t) + BT^'^{t) + V2ATQ'^{t), (2.4) 
tT„^'^(t) = AT„^;f(t) + BT„^'^(t) + AT„^_:f(t), n>2, 

where A is a symmetric and non-singular p x p matrix and i? is a symmetric p x p matrix, and 
T^'^(t) = Ip. Similarly, the Chebyshev polynomials of the second kind are defined by the recursion 

t[/„^'^(t) = A^U^f.it) + BU^^^it) + AUt'lit), n > 0, (2.5) 

with initial conditions U^{^{t) = 0, UQ'^(t) = Ip. Note that in the case p=l, A = 1,B = the 
polynomials T^'^(t) and U^'^{t) are proportional to the classical Chebyshev polynomials of the 
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first and second kind, that is 

in,x r- , in,x sin ((n + 1) arccos 4) 

T.i'O t = y2cos narcos- ; ^7^'° t = ^\ , ^ . 

" ^ ^ ^ 2^' " ^ ^ sm(arccos|) 

In the following discussion the p x p matrices A and B will depend on a real parameter u, 
i.e. A = A{u),B = B{u), and corresponding matrix measures of orthogonality are denoted by 
^A{u),B{u) (for the polynomials of the first kind) and Wa(u),b{u) (for the polynomials of the second 
kind). These measures will be normalized such that 



dXA{u),B{u){t) = J dWA(u),B{u){t) = Ip- 



If /3i(m) < . . . < Pp{u) are the eigenvalues of the matrix B{u), then it follows from Lemma 2.4 in 
Duran (1999) that the matrix 

KAiu)Mu){^) - - zIpy/'A-\u){B{u) - zQ'/^ (2.6) 

can be diagonalized except for finitely many z e C\[/5i(m), l3p{u)]. Throughout this paper the union 
of the set of these finitely many complex numbers and the set [/?!, {u),j3p{u)] will be denoted by 
A (note that the set A depends on the parameter u > although this is not reflected in our 
notation). In this case we have 

KAiu)Mu)iz) = U{z)D{z)U-\z), z e C\A, (2.7) 

where 

D{z) = diag (Af (")'^^")(z), . . . , A^(")'^(")(z)) 

denotes the diagonal matrix of eigenvalues of Ka{u),b{u) (^) and U {z) is a unitary pxp matrix [note 
that the matrices U{z) and D{z) depend on the parameter u, although this is not reflected by our 
notation]. It is shown in Duran, Lopez- Rodriguez and Saff (1999) that under the assumption of 
a positive definite matrix A{u) and a symmetric matrix B{u) the measure Xa{u),b{u) is absolute 
continuous with respect to the Lebesgue measure multiplied with the identity matrix and has 
density 

dXAiu),B(n)i^) = A-^/\u)U{x)f{x)U^{x)A-^/\u)dx (2.8) 
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where T(x) := diag(tn{x), . . . ,tpp{x)) denotes a diagonal matrix with elements 



1 



7r,/4- A 



A(u).B(u) 



f, if A 



A{u),B(u) 



x) e (-2,2) 



tii(x) = ( i = l,...,p. (2.9) 

0, else 
For the sake of a simple notation this density is also denoted by Xa{u),b{u){x). 

Theorem 2.1. Consider the sequence of matrix orthonormal polynomials defined by the three 
term recurrence relation Ii2. where for all £ G Nq and a given u > 



lim An-e,k = -4(ti), 



lim B. 



n-l,k 



(2.10) 
(2.11) 



with non-singular and symmetric matrices {A[u)\u > 0} and symmetric matrices {B{u)\u > 0}. 
// there exists a number M > such that 



OO CO f 

IJ IJ L G C| det Rn,k{z) = I C [-M, M] 

fc=ln=0 ' 



(2.12) 



then the empirical measure 6n,k defined by \2.'J\] converges weakly to a matrix measure which is 
absolute continuous with respect to the Lebesgue measure multiplied with the identity matrix. The 
density of the limiting distribution is given by 



1 

u 







/ tr 


-Xa{s),B{,s) 


Jo 





ds. 



(2.13) 



where Xa{s),b{s) denotes the density of the matrix measure corresponding to the matrix Chebyshev 
polynomials of the first kind. 

Proof. To be precise, let {Rn,k{x))n>o {k G N) denote the sequence of matrix orthonormal 
polynomials defined by the recursive relation (12.11) and denote by Xn+i,k,i, ■ ■ ■ , Xn+i,k,m the different 
roots of the (n + l)th polynomial Rn+i,k{x). At the end of the proof we will show the following 
auxiliary result, which generalizes a corresponding statement for the case p = 1 proved by Kuijlaars 
and van Assche (1999). 
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Lemma A.l. // all roots of the matrix orthogonal polynomial Rn+i,k{x) defined by Ii2.1\) are 
located in the interval [—M,M], then the inequality 

holds for all vectors v G C^, and all complex numbers z E C\[— M, M]. Moreover, we have 



(2.14) 



1 



(2.15) 



for all vectors v G C^, z E C with \z\ > M. 



We now normalize the ortho normal polynomials Rn.k{x) such that their leading coefficients are 
equal to the identity matrix, that is 



Eo^ki^) ■■= Ip ; Rn,ki^) ■= M,k ■ ■ ■ An,kRn,ki^), U > 1. 



(2.16) 



Then a straightforward calculation shows that Rjl{x)Rj^i i^{x) = Rjl{x)Aj^i^kRj+i k{x) (j > 0) 
and we obtain 

n— 1 

Rn,k{^) = l[Rj,l{^)A,^,,kR,+iA^), n>0. (2.17) 



j=0 



This yields 



— log|det^„^fc(z)| = — ^log|det (i?J(z)Aj+i,fei?j+i,fc(2))| 
np np 

= - I log Idet (Ai,^]+i^kRlsn]+l,k{^)Risn] fc(^) 



(2.18) 



ds. 



If Vn,k,ii^)^ • • • 5 Vn,k,piz) denote the eigenvalues of the p x p matrix Rn,kiz)Rnlik^)A^ki^, then it 
follows 



mm 



v^Rn,k{z)R^Lki^)A„l,^v 



< 



Vn,k,i[^) 



< max 



v^Rn,k{z)R;;kkiz)A~l,^,v 



for alH = 1, . . . ,p [see Horn and Johnson (1985), p. 181]. With these inequalities and Lemma A.l. 
we have 



\det{Rn,k{^)RnllA^)Kll 



In 



i=l 



'nn,k,i[Z} 



> min \'nn,k,iiz) 



i=l 



(2.19) 



whenever \z\ > M, and 



\det{R^4z)R-l,^,{z)A-l,^,)\ < (mjix|r/„,,,,(^)|)' (2.20) 



< ( max \v'^Rn,kiz)R^l-^^ki^)A^l-^^f^v\y < 



KvTv=i^ ' ' "+1''^^^ "+1''^ 7 - dist(2, [-M,M])P' 
From (12.191) and (I2.20p we therefore obtain the estimate 

|log |det {An+i,kRn+iAz)Rnl{z)) 1 1 < max {p |log (dist(^, [-M, M]))| ,p |log(2|;2|)|} , (2.21) 

whenever \z\ > M. Now the representation (I2.18P and Lebesgue's theorem yield for \z\ > M 

lim — log|det^„ ,t(z)| =- / log|det( lim v4[,„]+i,fci?M+i,fc(2;)i?r;^i ^(2;)) (2.22) 

^->u up p Jq -j^^u 



1 /•! 

- / log I det {^A(su),B(sn){.z)) \ ds, 
P Jo 

— / log I det ($A(s),B(s)(^)) I ds, 
Jo 



where 



^A{u),B{u)[Z) ■- / [Z.Z6) 

denotes the (matrix valued) Stieltjes transform of the matrix measure Wa{u),b{u) corresponding to 
the matrix Chebyshev polynomials of the second kind, and we have used Theorem 2.1 of Duran 
and Daneri-Vias (2001) for the second equality, which describes the "ratio asymptotics" of the 
matrices Rn+i,k{z) R^\{z). In fl2.22p the convergence is uniform on compact subsets of C \ F, 
where the set F is defined by 

m>0 j>m 

Aj = {x E M \ det Rnj,kj{x) = 0} denotes the set of roots of the matrix polynomial Rn^^kj{x), and 
{nj,kj)j>o is the sequence along which the convergence is considered. The Stieltjes transform of 
the matrix measure Wa{u),b(u) has been determined by Duran (1999) as 

$A(n),B(.)(^) = \A-\u)iBiu) - zl,f'^ (-/, - ^J, - 4K-J^),s(„)(2:)) {B{u) - zI,Yl'A-'\u) 
= A-\u){B{u) - zIj,Y/^U{z)T{z)U-\z){B{u) - zIpY/^A-\u), (2.24) 

where the matrix Ka{u),b(u) is defined in (12.61) . T{z) := diag (tii(z), . . . , tpp{z)) is a diagonal matrix 
with elements given by 



tii{.z) 
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1, . . . ,p) and z E C \ {supp {Wa{u),b{u)) U A}. Here and throughout this paper we take the 



square root y/w such that \w — — 4| < 2 for w G C \ [—2, 2], and consequently the function 
w — y^w"^ — 4: is analytic on C \ [—2,2]. Note that for z G C \ supp (Wa{u),b{u)) it follows that 
^ 2. Observing (12 .24^ this implies for the Stieltjes transform 

log|det($^(„),BH(^))| = \og\det{A-\u)T{z)iB{u) - zlp)\ (2.25) 



log I det{A'~\u)T{z)D{z) \ = logfQ 



ctj (m) 



^log^ V 



2aj(u 



where «!(«), . . . , ap(u) denote the eigenvalues of the matrix A{u). Now define 

_^Aiu)Mu) _ ^(A^^W'^H(^))2_4 
2a,- (m) 

then it follows from Kato (1976), p. 64, that the eigenvalues X'^^^^'^^^\z) are holomorphic functions 
on C \ A and we obtain for any 

z E Gq := C\ {supp {Wa{u),b{u) 



that 



^ 1 

+ 1 

^ 2/j(2;)a,(M) V ^(;Af"^'^^"^(;z))2 -4 

= y , ^ ( ^Af")'^(^)(z) +4Af")'^(")(;.)^ = 

^ 2J(Af")'^(")(z))2-4 ^ / 



where the function is defined by 

_^Aiu)Mu)^^^ _ ^(a;iw.^(«)(^))2_4 

_ . . , 7 = 1 



2aj[u) 



fjiz) = 7^ — —\ ' j = 1, • ■ ■ 

10 



This implies that the function f{z) is holomorphic on Gq, and it follows for z G Go that 



= y ^ ^^""^ -^^ (2.26) 

In the following discussion define 

W'iz) := J log^l^tr [XAiu)MuM dt (2.27) 

as the logarithmic potential of the measure whose density with respect to the Lebesgue measure 
is given by tT[XA{u),B{u)if)\- Then U^{z) is harmonic on 

Gi := C \ SUpp {tT[XA{u),B{u)]) 

[see e.g. Saff and Totik (1997)]. In the following we show that the function / satisfies Re / = 
+ Ci, where Ci G C is a constant. For this purpose we note that the function 

,i.):=±uH.)-,£u-i.) (2.28) 

is holomorphic on Gi (note that g satisfies the Cauchy-Riemann differential equations because the 
logarithmic potential is harmonic) and satisfies for all z E Gi 



9{^) = - I ^3^tr[X^(„),B(„)(t)]rft (2.29) 

tr \XA{u),B(u){t)\ 



where Ga(u),b{u) denotes the Stieltjes transform of the matrix measure corresponding to the density 
XA{u),B{u){t)- In order to find a representation for the right hand side we note that the function 
Ka(u),b{u){z) defined in (12. 6p is holomorphic on C \A, and observing the representation (12. 7p we 
obtain for all z G C\A 

-^{B{u) - zI,)-'/'A-\u){B{u) - zl,y/' - ^{B{u) - zIj,y/'A-\u){B{u) - zQ-'/' 
= K'a^uIbU^) = U'{z)D{z)U-\z) + U{z)D'{z)U'\z) + U{z)D{z){U~\z)y . 
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the identities 



A-\u){B{u) - zl^y/^ = {B{u) - zIp)-'/^U{z)D{z)U-\z), 
{B{u) - zIpY^^A-'iu) = U{z)D{z)U-\z){B{u) - zI^Y^'^ 



and {U~\z)yU{z) = -U-\z)U'{z) it follows that 

-^U-\z){B{u) - zQ-^U{z)D{z) - ]^D{z)U-\z){B{u) - zQ-^U{z) (2.30) 
= U-\z)U\z)D{z) + D\z) - D{z)U~\z)U'{z), 

which yields for the diagonal elements of the matrix U^^{z){B{u) — zIp)^^U{z) 



[U-\zmu) - zl.ruiz)],, = (-[^-^(.)^-(.)kAf")-^(")(.) - ^Af ")'^(")(.) 



Af")'^(")(^)[f/-^(.)f/'(.)],, 



for all z G C\A, j = 1, . . . ,p. From Zygmunt (2002) and (12.241) we have for the Stieltjes transforms 
^A{u),B{u) and Ga{u),b{u) corresponding to the matrix measures Wa{u),b{u) and Xa{u),b{u) the matrix 
continued fraction expansion 



GA(u),Biu){z) = {zIp-B{u)-V2A{u){zI -B{u)-A{u){zlp-B{u)-... 
... A{u)y' V2Aiu)y' 
= {zip - B{u) - 2A{u)<^A{u),Biu)iz)Mu)y^ 
= {B{u) - zIp)-^/^U{z)f{z)U-\z){B{u) - zlp)-y\ 

where z G G, the set G is given by 

G := C\{supp {Wa(u),b{u)) U supp {Xa(u),b(u)) U A}, 

and the diagonal matrix T{z) is defined by 

.A{u),B{n)( \ .A{u),B{u),. 

f{z) := diag(tn(^), • • ■ = diag( ^ > , , " ^ ' =). 

\^^^Aiu),Biu)^^^y _ 4 ^^^Aiu)Mu)^^^y _ 
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This yields 



g{z) = -tv[GAiu)MuM = -J2[U-\z){B{u)-zI,)-'U{z)y^ 

i=i 



p ■.A(«),_B(«)/ N 

= \^ ^ lJ. (2 31) 

for all z & G. Consequently, we have for all 2; G G 

^ I ^yMumu) ( ,^^2 _ A dx dy 

For /i := Re / it follows that 



,=1 W(A;^"^'^^«^(z))2-4 



da; ay ax ay 
so that j^ih - U^) = and ^{h - U^) = 0, which implies for all 2; G G the identity Re f{z) = 
U^{z) + ci for some constant Ci G C. Therefore it follows that 

Re (/(z)) = Re (log(det($AH,BM(^)))) = log | det(<l>A(„),B(n)(^))| = U"" (z) + Ci (2.32) 

for any z & G. Observing fl2.22p we finally obtain for all 2; G G 

1 1 /" 1 

li"^ -l°g |det^„,,(^)| = -- / / log^^t^ [^A(s),B(.)(t)] cit - c (2.33) 



1 1 r 

log^ 7\- / tr 



-^A(s),B(.)(t) 



dt — c 



= -U^iz)-c, 

where If^ denotes the logarithmic potential of the measure a with Lebesgue density defined 
in (12.131) and c G C is a constant. Let 

denote the logarithmic potential of the measure 5n,k defined in (12. 3p . From (I2.33P we obtain for 
all z G G 

1 1 1 
lim U^"'''(z) = hm — > log-; -r- = lim — lo£ 



\z-Xn,k,j\^' TT^^np \det{Rn,k{z))\ 

f/"(z) + c. (2.34) 
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The measures in the sequence ((5„^.^fcJjgN have compact support in [— M, M], consequently, (5„^ ,A:j )jeN 
contains a subsequence which converges weakly to a limit fi with supp (/i) C [— M, M]. Therefore 
we obtain from fl2.34p 

U^'iz) = U^z) + c, zeG, \z\ > M, 

and the assertion of the theorem follows from the fact that the logarithmic potentials are unique 
[see Saff and Totik (1997)]. 

Proof of Lemma A.l. Let Xn+i,k,i, ■ ■ ■ , Xn+i,k,m denote the different roots of the matrix orthog- 
onal polynomial Rn+i^k{x) with multiplicities ii, . . . ,im- Then we obtain from Duran (1996), p. 
1184, the representation 

ILA^)R-lJz)A-l,, = J2 ""^^-^ (2.35) 
where the weights are given by 

Cn+l,fc,j^„+l^fc = Rn,kiXn+l^k,j)^n+l,k,jRn,ki-^n+l,k,j) 1 (2.36) 

^k 

Tn+i.fcj = — — , ,.,(1) , (Adj(i?„+i,fc(t)))^^^"-^^(x„+i,fcj)Q„+i,fc(a;n+i,fcj) , 

(det(i?„+i,fc(t)))^'^^ (x„+i,fcj) 

and Adj(A) denotes the adjoint of the p x p matrix A, that is 

A Adj(A) = Adi{A)A = det{A)Ip . 

In (12.371) the matrix polynomial 

Qn,k[x) = / 7 aSfc(t) 

J X - t 

denotes the first associated matrix orthogonal polynomial and the matrices Tn,k are nonnegative 
definite and have rank [see Duran (1996), Theorem 3.1 (2)]. From Duran and Daneri-Vias 
(2001) we have 

m m 

^ C'n+l,fcj^„+l,fc = '^^Rn,k{Xn+l,k,j)^n+l,k,jRi,k{^n+l,k,j) (2.37) 
3=1 3=1 

Rn,ki^)^kit)Rn,ki^) = Ip- 

For any z G C\[— M, M] we obtain the estimate \z — Xn+i,k,j\ > dist(2;, [— M, M]) for all j = 
1, . . . ,m, because all roots of the matrix polynomial Rn+i,k{z) are located in the interval [— M, M]. 
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Note that by the representation (I2.36P the matrix Cn+i,k^nl-i k nonnegative definite and the 
representations (12.351) and fl2.37p yield for any |2;| > M and v & 

\ Tr> I \ D-1 ( \ I T ri+l,fc 

m T 



< 



1^ "^n+lj/cjl 



dist(;2, [-M,M]) t.„+i,fcjA„+i,fcT; = dist(;2, [-M, M])' 



which proves the first inequahty of the Lemma. For a proof of the second part we note that for 
1^1 > M we have |^i:i±M:i| < \ for all j = 1, . . . , m, which gives 



Re r . ■ > - 



for all j = 1, . . . , m. With this inequality we obtain 



I I j=l Z 



m T 



1 t;" G, 



m 



which proves the second inequality of the Lemma. □ 



j=l 



3 Strong and weak asymptotics for eigenvalues of random 
block matrices 

We now consider the random block matrix defined in equation flL3l) of the introduction and 
denote A["'^'' < •■■ < Xi?'^'' as its (random) eigenvalues, where n = mp; m,p E N. Similarly, 
define x^^'^^ < ■ ■ ■ < (deterministic) eigenvalues of the matrix orthonormal polynomial 

R';^]nix) E which is defined recursively by R'^lix) = Ip, ^i^J,„(x) = 0, 
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(m > 0), where the p x p block matrices A[^^ and B^^^ are given by 

f l)P + l)7p l)P + 2)7p-i l)P + 3)7p-2 

l)p + 2)7p-i l)P + 2)7p V'((«-l)P + 3)7p-i 



/ip7l \ 



V" («p - 1)72 



^{ip - l)7p- 



V^PTp / 



V*P7p- 



(3.2) 



i > 1, and 



(p) 



1 

71 



v' (jp + 1)71 



V (ip + 1)71 V (*P + 1)72 
V(»P + 2)71 



V(^P + l)7p-l \ 
V(*P + 2)7p-2 



V(«P + l)7p-2 
V V(«P + l)7p-i 



^((i + 1)P - 2)71 ^({i + 1)P - 2)71 

+ 1)P - 2)72 ^/({i + 1)P - 1)71 



(3.3) 



i > 0, respectively, and 71, . . . ,7p > are given constants. Our first result provides a strong 
uniform approximation of the ordered random eigenvalues of the matrices G^'' by the ordered 
deterministic roots of the matrix orthogonal polynomials Rm]n{x). 

Theorem 3.1. Let A^"'^^ < ■ ■ • < A^"'^^ denote the eigenvalues of the random matrix Gn^ defined 
by l{1.3\) and x^"'^'' < ■ ■ ■ < Xn"'^'' denote the roots of the matrix orthonormal polynomial Rm]n{x) 
defined by the recurrence relation liS. then 



max 

i<i<" 



^(n,p) _ ^(n,p) ^ [log^]l/2^ Vn>2, 

where S is a random variable such that S < 00 a.s. 



(3.4) 



Proof. For the proof we establish the bound 



P i max 

l<j<n 



\("'P) _ ™("'P) 



> e > < 2n{p + 1) exp 



(3.5) 



then the assertion follows along the lines of the proof of Theorem 2.2 in Dette and Imhof (2007). 
First we note that the roots of the mth matrix orthonormal polynomial Rm]n{x) are the eigenvalues 
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of the tridiagonal block matrix 



^l,n Bin ^2,n 



\ 



■2,n 



Br. 



■2,n 



V 



(3.6) 



V V / 



where the blocks v4.^, i = 1, . . . , n/p 



1, and 5^2, i 



0, . . . , n/p — 1, are defined by (I3.2p and 



(I3.3p . respectively. Moreover, by interchanging rows and columns of the matrix Fn'' it is easy to 
see that the matrix 



(P) 



(P) 
l,n 



2,n 



n(?') 7?(p) n(p) 

^^-2,n ^^-2,n -^^-l,n 



-l,n 



E 



(p) 



-l,n 



(3.7) 



with symmetric blocks 



/ 

V71 (" - ip - 1) 



V2 



1/71 (n - ip - 1) v'72(n - ip - 2) 
^7i{n - jp - 2) 



V7p-2(n- (i + l)p + 2) 
V V7p-l(n- (i + l)P+l) 



^ - 1, and 



V7p-i("- (i + l)P + l) ^ 
^7p_2(n- (i + l)p + l) 



^71 (n- (j + l)p + l) ^-yi{n - (i + l)p + 1) 

V72{n- (j + l)p + l) V7i(™- (i + l)P+l) 



^7p_i(n - ip) V'7p(n-«P- 1) V'7p-i{" - «P- 1) 



(p) 



1 

71 



^71 (n - ip) 
^72(n- ip - 1) 



^72 (n - ip) 
V V7i (" - «P) 



V7p-l{n- (j + l)p + 3) v'7p("- (i + l)P + 2) v'7p-i("- (i + l)p + 2) 



V7p-2(n- (i + l)p + 3) V7p-i("- (i + l)P + 2) V7p ("-(« + 1)P + 1) / 
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(z = 1, . . . , ^ — 1) also has eigenvalues x^"'^^ < ■ ■ ■ < Xn'''^^ ■ From Horn and Johnson (1985) we 
therefore obtain the estimate 



max 

l<j<n 



n n 



where for a matrix A = (aij)i ,=i 



_ = max > a,-,: 

l<i<n — 



denotes the maximum row sum norm. If 



Z„ := max { max I A'', ! , max —= 
_i<i<" i<i<n-i V2 



max 



l<j<n-p ^ 



i7p y Jlp 



then it follows for i = 1, . . . , p, n — p + 1, . . . , n by a straightforward but tedious calculation 



and for i = p + 1, . . . , n — p 



which yields observing (13. 8p 



and 



<3pZ., 



max 

l<i<n 



P { max 



3p 



From Dette and Imhof (2007) we obtain the estimates 
P 



> 



y/2 SpJ ll<j<n 

and a combination of these inequalities with (13. 9p yields 



< 2e-'"/^P" , P\ max |Ar,| > < 2ne-''/^^P' 



e 

3p- 



P < max I A,- — X 



I l<i<n 

which completes the proof of the theorem. 



>£|<p|z„>^|< 2n{p + l)e-^'/i^P' 



(3.9) 



□ 
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Theorem 3.2. Let X["''^^ < • ■ ■ < An"'^'' denote the eigenvalues of the random matrix -^Gn\ 



where Gn^ is defined in U.3\) and 71, . . . , 7p > are chosen such that all blocks Df^ in the matrix 



(p) 



Fn'' defined by 1^3. ?| j are non-singular. If 

1 " 



(3.10) 



denotes the empirical distribution of the eigenvalues of the matrix -^Gn \ then cr„ converges weakly 
to a measure, which is absolute continuous with respect to the Lebesgue measure. The density of 
this measure is given by 



fit) 



tr [Xj^ip 



)(s),B(p)(s) 



(t)] ds, 



(3.11) 



where -^a(p)(s),s{p)(s)(^) denotes the Lebesgue density of the matrix measure corresponding to the 
Chebyshev polynomials of the first kind defined in \2.4\) with matrices 



( 




7p VV^ 



7i 
7^ 




72 

tT 



71 



71 




VtT 
V72 Vti 



3Xp 



(3.12) 



Proof. Let i?}^„ (a;) (m = n/p) denote the orthonormal polynomials satisfying the recurrence 
relation with coefficients A[^1 = ij^^ and 5^^^ = ^ where ij^^ and B^^} are defined 



by fl3.2p and (13. 3p . respectively, then we have 



X 



(n,p) 



n 



where the matrix polynomial Rm]n{x) is defined by (13.11) with corresponding roots x^"''^'', and x^' 



„(".p) 



denotes the jth root of the matrix polynomial Rm]n{x). From the definition of the matrices AfJ^ 



and in (D and ([SJD we have 



ip) 



lim A_ 

n— »oo p 



-l,n 



1 

7! 



7i 

72 



72 



71 



( 



hm Elf) 1 „ = ^ 



71 




v/7^ VTp^ 



72 

71 



tT 

72 Jli 



7i 




J 



and Gerschgorin's disc theorem implies that all roots of the polynomials R^m,n{x) are located in a 
compact interval, say [— M, M] [see also the proof of Lemma 2.1 in Duran (1999)]. Moreover, we 
have for any £ G No 

lim SH„ = 



upA^P'^ =: A'^P\u), 



n 

where m > and the matrix A^'p\u) is non-singular by assumption. Consequently, Theorem 2.1 
is applicable and yields (note that lim„^oo('^/p)/'T- = that the empirical distribution of the 
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roots of the matrix polynomials -Rm)n(a;) (m = n/p) 

1 " 

= -Y^^x'^-'^^ (3.13) 



n — ' -^j 



converges weakly to the measure with Lebesgue density / defined in fl3.1ip . Next we use Theo- 
rem 3.1 which shows that 



max I A^^) - xf'^ \ = ^ max | A^^) - | = O ( f '^') (3.14) 



almost surely. Therefore we obtain for the Levy distance L between the distribution functions F^r^ 
and Fs^ of the measures (T„ and 5„ 

L'iF^., FsJ < - E I - P = O (^) (3.15) 

almost surely [for the inequality see Bai (1999), p. 615]. Consequently, it follows that the spectral 
measure cr„ of the matrix -^Gn^ converges also weakly with the same limit as 6n, that is the 
measure with Lebesgue density / defined by 03. lip . 

□ 



4 Examples 

We conclude this paper with a discussion of a few examples. First note that in the case p = 1 
Theorem 3.2 yields for the limiting distribution of eigenvalues of the matrix -^Gn'' the Wigner's 
semicircle law with density 

which have been considered by numerous authors. 

Next we concentrate on the case p = 2, for which it is easily seen that the matrix D^"^^ in 03.81) is 
non-singular whenever 72 7^ 7i. In this case the density of the limit of the spectral measure is a 
mixture of two arcsine densities and given by 
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Figure 1: Simulated and limiting spectral density of the random block matrix Gn j \fri in the case 
p = 2, 7i = 2, 72 = 8. In the simulation the eigenvalue distribution of a 5000 x 5000 matrix was 
calculated (i.e. m = n/p = 2500). 

Figure 2: Simulated and limiting spectral density of the random block matrix G^n^ j ^Jn in the case 
JO = 2, 7i = 1, 72 = 100. In the simulation the eigenvalue distribution of a 5000 x 5000 matrix 
was calculated (i.e. m = n/p = 2500). 

where 

In Figure 1 and 2 we display the hmiting spectral density corresponding to the case 71 = 2, 72 = 8 
and 7i = 1,72 = 100, respectively. The left part of the figures shows a simulated histogram of 

1 (2) 

the eigenvalues of the matrix Gn for n = 5000, while the right part of the figures shows the 
corresponding limiting distribution obtained from Theorem 3.2. 

If p > 3 the general formulas for the density of the limit distribution are too complicated to be 
displayed here, but it can be shown that 

,=1 vry/4 - ^xf Hu),BM(u)^^-^y {-2<A, {t)<2} 

if the matrix A^p\u) is positive definite. This identity follows in a similiar way as fl2.3ip . In Fig- 
ure 3, 4 and 5 we show a simulated histogram of the eigenvalues of Gn^/^Jn and the corresponding 
density of the limit distribution obtained from Theorem 3.2 in the case 71 = 72 = 4, 73 = 100; 
7i = 1, 72 = 4,73 = 25 and 71 = 1, 72 = 100,73 = 200, respectively. 
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der forschungsbereich Tr/12, Fluctations and universality of invariant random matrix ensembles 
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Figure 3: Simulated and limiting spectral density of the random block matrix Gn / \fn in the case 
p = 3, 7i = 72 = 4, 73 = 100. In the simulation the eigenvalue distribution of a 5001 x 5001 
matrix was calculated, (i.e. m — n/p — 1667) 

Figure 4: Simulated and limiting spectral density of the random block matrix C't^ j ^Jn in the case 
p = 3, 7i = 1, 72 = 4, 73 = 25. In the simulation the eigenvalue distribution of a 5001 x 5001 
matrix was calculated (i.e. m = n/p = 1667). 
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